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Abstract 

The effect of dynamical short-range correlations on the general- 
ized momentum distribution n{p, Q) in the case of Z = N, ^-closed 
shell nuclei is investigated by introducing Jastrow-type correlations 
in the harmonic-oscillator model. First, a low order approximation is 
considered and applied to the nucleus ^He. Compact analytical ex- 
pressions are derived and numerical results are presented and the effect 
of center-of-mass corrections is estimated. Next, an approximation is 
proposed for n{p, Q) of heavier nuclei, that uses the above correlated 
n{p, Q) of ^He. Results are presented for the nucleus ^^O. It is found 
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that the effect of short-range correlations is significant for rather large 
values of the momenta p and/or Q and should be included, along with 
center of mass corrections for light nuclei, in a reliable evaluation of 
n{p,Q) in the whole domain of p and Q. 

Keywords: Generalized momentum distribution; Density matrices; Quasi- 
elastic scattering; Final state interactions; Momentum distributions; Short- 
range correlations. 
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1 Introduction 

In a system of A {A > 2) identical particles described by a unit-normalized 
state the generalized momentum distribution n{p,Q) (GMD) is defined 

by EEIE] 

k,s,s',t,t' 

By introducing the density fluctuation operator = J^j: s t ^\+q s t^k s t 
and the single-particle momentum distribution n{p), definition ([T]) is recast 
as 

n{p.Q) =Mmp^Y.^l_^,s,^,^,s,m-n{p)] . (2) 

s,t 

The role of GMD in final state interactions becomes evident from this ex- 
pression, since the first term on the right is a transition matrix element for 
scattering a particle out of the orbital (p, s, t) into the orbital {p—Q, s, t), this 
process being introduced by a spin-isospin independent density fluctuation 
of wave vector Q. 

The quantity n{p, Q) is connected to the two-body density matrix (2DM) 
in momentum space n2{pi, P2] Pv , P2') through the relation 

n{p,Q) = J n2ip,k;p-Q,k + Q)d^k . (3) 

Introducing the two-body density matrix in the coordinate space and its half 
diagonal version P2h{'^i,'^2]'^r) and performing a Fourier transform in the 
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coordinates fl — ry and ry — we obtain for GMD 




We shall take ^ = c = 1, so that the momentum has the dimension of inverse 
length and n{p) and the GMD have the dimension of (length)'^. The normal- 
izations adopted here are such that / n2{'Pi,P2]Pi-,P2)^Pi^P2 = — 1) 
and / n{p)(^p = A. 

Regarding the above definitions, in Refs.[l-3] no distinction was made 
between "laboratory" momenta and intrinsic momenta, i.e., momenta with 
respect to the Center of Mass (CM) frame. For infinite systems, such as 
considered in Refs.[l,2], CM correlations are not relevant. For large self- 
bound systems like heavy nuclei, they are also not very important. For 
small systems, however, like the nucleus ^He, one has to be careful with the 
definitions and relations given above. In particular, the momenta in Eqs. 
([T])-(jl]) should be interpreted as momenta with respect to the center of mass 
of the system. We shall return to this point later. 

The generalized momentum distribution n{p, Q) has some important for- 
mal properties that result from the corresponding properties of the 2DM in 
momentum or coordinate space. In particular, the sequential relation be- 
tween the half-diagonal 2DM and the one-body density matrix (IDM) yields 
the relation between the GMD and n{p), namely 

n{p,Q = 0) = {A-l)n{p} . (5) 

— * 

In addition time-reversal invariance implies that n{p, Q) is symmetric with 
respect to the variables p and w {w = p — Q) namely 

n{p, Q) = n{w, -Q) . (6) 

The quantities n{p,Q) and P2hifi,r2]rii) are important descriptors of 
nucleon-nucleon correlations in the nuclear medium representing the next 
stage of complexity beyond the single-particle momentum distribution and 
the one-body density matrix. Therefore, the last years interest has risen for 
their study. One main reason is the need to properly analyze recent and 
future experiments of inclusive character, such as (e,e'), (p,p') reactions, as 
well as of exclusive character such as (e,e'N), (p,2p), (7,N) and (e,e'2N), 
(7,2N) and to extract reliable values for the momentum distribution, the 
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one- and two-particle spectral functions, the transparency and other quanti- 
ties To achieve these goals, one must take into account the final-state 
interactions (FSI) of the struck nucleons as they propagate through the nu- 
clear medium. The half-diagonal 2DM appears in almost all quantitative 
microscopic post-mean field treatments of the FSI (see for instance, in the 
case of inclusive (e,e') scattering [21] and references therein). The GMD, as 
shown before, (Eq. ([2])), is directly involved in FSI mechanisms. Finally, in 
analogy with other quantum many-body systems [251 ISSl [27], the functions 
P2h{^i,r2',ri') and n{p,Q) are expected to enter in fundamental sum rules 
that furnish insight into the nature of elementary excitations of the nuclear 
system. 

Initially, both the half-diagonal 2DM and the GMD were investigated 
within the context of final-state effects in inelastic neutron scattering from 
quantum fluids (i.e. liquid He). With regard to nuclear systems, calculations 
based on the work of Ristig and Clark [T] within the context of variational 
theory have been performed for Jastrow-correlated inflnite nuclear matter. 
This has been done for the GMD using low-order cluster truncations [28] 
and later for the GMD and the half-diagonal 2DM using Fermi-hypernetted- 
chain procedures [21 12H] • These results clearly show the effect of short-range 
correlations (SRC) in the nuclear medium and could be used by means of a 
suitable local density approximation for the approximate evaluation of the 
GMD of flnite nuclei. At the same time, the development of rather simple 
expressions for the latter which could be easy to use, led to the study of 
GMD of closed-shell nuclei within the independent-particle model and to the 
extraction of closed analytical expressions using a harmonic oscillator (HO) 
basis [3]. Only the proton contribution to GMD was considered and the 
approach applied was an extension of the one applied in Refs. [30l [311 [32] 
for the calculation of the charge form factor, the nuclear charge, matter 
and momentum distribution and the one-body density matrix. Recently, the 
same approach has been used for the calculation of the two-body momentum 
distribution [33j. The above results for GMD of nuclei exhibit interesting 
features rising from the flnite size and the Fermi statistics and are expected to 
be valid in certain regions of momenta p and Q where dynamical correlations 
do not play a signiflcant role. 

In this paper, we extend the above study of GMD in Z = N, £-closed 
shell nuclei by considering the proton and neutron contribution and including 
Jastrow-type correlations via the lowest term of a cluster expansion. This so- 
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called low-order approximation (LOA) has been exploited for the IDM 
and 2DM by Bohigas and Stringari [33] and Dal Ri, Stringari and Bohigas 
[36j and has been used to study single-particle nuclear properties |35j-[39] 
as well as two-body nuclear properties [33l HHl HI]. Recently it has been 
extended for the case of realistic interactions and corresponding correlations 
and applied in the calculations of the ground-state energies, densities and 
momentum distributions of ^^O and ■*°Ca [121113]. In our application of LOA 
to the evaluation of GMD we have used central single-Gaussian correlation 
functions, /(r) = 1 — ce~^ , and we have performed explicitly the calcula- 
tion of GMD for "^He. Due to its high central density (almost 8 times nuclear 
matter density) the nucleus ^He is a particular appropriate system to search 
for SRC. We find significant deviations from the independent-particle model 
picture due to SRC for rather large values of momenta p and (or) Q. We 
also examine the effect of Center of Mass corrections and we find that it is 
quite significant and should be taken into account. The calculation for heav- 
ier ^-closed nuclei is more complicated. We have developed an approximate 
scheme which makes use of the above calculated GMD of "^He and which is 
likely to be valid for not very large mass number A. Results are extracted for 
^^O. It should be emphasized that the approach presented for the calculation 
of the GMD using Jastrow-type correlations is simply an exploratory one, 
aimed at guiding realistic calculations which include a full correlation oper- 
ator. A similar procedure has been followed in the study of other quantities 
including the momentum distribution. 

In section 2 a brief outline of the calculation of the GMD for Z = N 
^-closed shell nuclei in the independent-particle shell model with harmonic 
oscillator wave functions is given. In section 3 a short presentation is made of 
our estimation of CM corrections in the evaluation of GMD. In section 4 the 
effect of SRC on the GMD is explored by including Jastrow-type correlations 
via LOA and results for the case of "^He are extracted. Subsequently, the 
approximation for the evaluation of GMD of heavier nuclei that uses the 
GMD of ^He is presented and it is used for the evaluation of GMD of ^^O. 
Finally, in section 5 a summary, conclusions and hints for further development 
are given. 
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2 GMD of ^-closed nuclei in the harmonic os- 
cillator model 

We consider finite nuclei in their ground state within the independent-particle 
model. Using the relation of GMD to the 2DM in the momentum space 
^2{Pi,'P2]Pi',P2') (Eq- (J3])) and the corresponding expression of the latter for 
a system of noninteracting fermions we derive the following expression 

n{p, Q) = AF{Q^)nMp- Q) + n'^'ip, Q) , (7) 

where F{Q'^) is the elastic form factor, ni{p,p') is the IDM in momentum 
space and rf^{p,Q) is the exchange term, arising from the statistical corre- 
lations among the noninteracting fermions generated by the Pauli exclusion 
principle. In the case of spin- and isospin-independent Hamiltonians rf^{p, Q) 
equals 

n'\p,Q) = -\j niip,k + Q)n^{k,p-Q)d^k , (8) 

{u the degeneracy due to spin-isospin). The evaluation of the proton contri- 
bution to GMD by means of Eqs. ([7]), ([8]) using harmonic-oscillator single- 
particle states - but ignoring the Coulomb potential among the protons, the 
effects of the center of mass motion and the finite size of the nucleons - has 
been presented in detail in Refs. [21 E2] in the case of j-closed shell nuclei. 
Here, we present the expressions in the case of ^-closed shell nuclei, since 
we will refer explicitly to the nuclei ^He and ^^O. Numerical results will be 
presented for Q parallel to p {Q = QpP, w = p — Q = {p — Qp)p = Wpp) and 
for Q perpendicular to p . We have for the elastic form factor F{Q'^) pi] 

^ A=0 

where b is the harmonic oscillator parameter, A^max = (2?T- + ^)max is the num- 
ber of energy quanta of the highest occupied proton level, and the coefficients 
dx are rational numbers varying with Z. Their expressions as well as their 
values for the ten lowest ni levels are given in App.A of Ref. [3]. The IDM 
in the case of p' parallel to p {p' = p'^p) is given by 

?j3 -/Vmax -^max 

Mp^p'p) = -ij-,e-^'''/'e-^'^'''/' J: {phY E W'K,,' • (10) 
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The coefficients i^^^/ are discussed in the App. B of Ref. From expres- 
sion (II up the expression of the spherically symmetric nucleon momentum 
distribution n{p) is obtained [30] 

n{p) = n{p) = ^e-^^^^ ^ (p6)^^2/, . (11) 

^ A=0 

Inserting expression (fTOj) in Eq. ([8]) the corresponding expression for the 
exchange term rf^ is derived 

ax 

ax 

5: {w,hY' y: iQpbyc,,,, . (12) 

li'=0 p=0 

The coefficients C^/^'p are equal to zero if /i + /x' + p =odd and are discussed in 
App.C of Ref. Using the values listed there one can calculate the GMD 
of ^-closed shell nuclei up to Z or equal to 40. Expressions for rii and n^^ 
for Q not parallel to p can be found in Refs. [3l [32] . 

The GMD calculated via the above expressions in the case of ^He, ^^O 
and ^°Ca, as a function of Qp for given p, exhibits a bump centered at Qp = 
for p = and shifted to higher values of Qp for p > 0. A negative part in 
the GMD appears at positive Qp, for nuclei heavier than ^He, arising mainly 
from the term rf^{p,Qp). A comparison with respective results for an ideal 
Fermi gas and for nuclear matter shows that the positive bump and negative 
part are bulk properties of the GMD due to Fermi statistics [3j. 



3 Center of Mass Corrections 

As mentioned in the Introduction, special care should be taken when dealing 
with finite self-bound systems like nuclei, especially small ones like '^He. The 
CM motion in such cases cannot be ignored. The wave functions which are 
used in the independent particle model (but even in theories which take also 
dynamical correlations into account e.g. Brueckner-Hartree Fock, Variational 
Monte Carlo) satisfy the Pauli principle but not the translation invariance. 
As a consequence, they contain spurious components which result from the 
motion of the CM in a non free state. Effects from these (also know as CM 
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correlations) are found in the calculation of almost every observable and make 
impossible to extract information for the intrinsic properties of nuclei directly 
from experimental data. We will follow Ref. |S| for the evaluation of CM 
effects on the GMD (a brief history of the CM problem is found there). In 
Ref. IH] the two-body density matrix and the two-body momentum distri- 
bution have been studied by using the fixed CM approximation to construct 
the intrinsic wave functions, the Jacobi variables to define the correspond- 
ing intrinsic operators and an algebraic technique based upon the Cartesian 
representation of the coordinate and momentum operators to evaluate the 
expectation values involved. We proceed accordingly for the evaluation of 
the intrinsic GMD. 

The momenta in Eqs. should be interpreted as momenta with 

respect to the center of mass of the system. The two-body density matrix, 
such as that used in Eq. ([3]), should describe the transition matrix element 
of two nucleons from intrinsic momenta 

PA - P/A = p, pA-i - P/A = k 

to respective intrinsic momenta p', k' . With pi we denote the momentum 
of the i—ih nucleon with respect to the artificial mean-field center and with 
P = J2t=i'Pi the CM momentum. Following the notation of Ref. [H] we 
write the relevant operator in the form (hats denote operators) 

n2{p,k;p',k') = A{A-1)\pa- P/A = P)(pa- P/A = p'\ 

®\pa-i - P/A = k)(pA-i - P/A = k'\. (13) 
Switching to Jacobi momenta, as defined in Ref. [44J, we have 

n2{p,k;p',k') = A{A - l)\f]A_^ = p){f]^_^ = p'\ 

^fvA-2 = k + j^p){ffA-2 = k' + j^p'\. (14) 

In Ref. j44] a different intrinsic TBMD operator n^^l was defined formally in 
terms of Jacobi momenta, with a somewhat different physical interpretation. 
The two operators, n2 and n^'^\ are connected via the relation 

n2(p, k; p', k') = n'^l (p, k + -r^p; p', k' + -j^p')- (15) 
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Notice that both operators depend only on intrinsic (Jacobi) variables. Since 
in the simple HO model the CM wave function factorizes into an intrinsic 
and a CM part, there is no need to project out the CM component of the 
wave function (by means of the fixed-CM or other method). 

The matrix elements of the operator 77,2 can be calculated using the for- 
malism of Ref. For ^He, i.e., in the state |0s^), its expectation value is 
given by 

n2(p, k] p' , k') = [ e -4-2 2 % a-2(p^+pI')\ (16) 

TT-^ \A — 2/ 

From 71,2 (p, k;p', k') we can evaluate the generalized momentum distribu- 
tion n(p, Q) by setting p' = p — Q and k' = k + Q and integrating over /c, 
according to Eq.(3). We find 

n{p,Q) = ^ (-A_)''\-^M-fQ)^\-¥^Q^^\ (17) 



7r3/2 \A - 1 

which can be written also in terms of w = p — Q as 



<P^Q) = % (-^)'''e-^(-^^-^)^^e-i^'3^^l (18) 



The corresponding expression for the intrinsic momentum distribution 
n{p) is [B] 

The generalized momentum distribution and the momentum distribution of 
^He in the HO model without CM corrections denoted by n^{p,Q]b) and 
n^{p;b) respectively are given by (see Eqs.(12), (11)) 

n'ip, Q; b) = l|!e-^^''V2,-^^.V2,-Q^^.V4 (20) 



n'm = ^,e-^^'^- (21) 



4b 



3 



vr 



3/2 



We have written in Eqs. f|T6|) - f|T9|) the different coefficients in terms of A and 
not 4 to point out a trend in A-dependence of the effect of CM motion. 

The variables p'^ and appear to scale due to the CM correlations by 
the factor (1 — \)~^, i-e. the inverse Tassie-Barker factor (TBF), just like the 
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variable in the momentum distribution [HI US] . The variable scales 
by the factor (1 — -^^)- These observations can provide an approximate 
scheme for introducing CM corrections to the GMD of ^He also when SRC 
corrections are considered, as we discuss in the next section. 

In heavier nuclei, CM corrections should be less important. To illustrate 
this, we note that the GMD in the simple harmonic-oscillator model for any 
nucleus is given by the same exponential as above, times a polynomial - 
see Eq. ( !T2|) and Ref. [3]. After including CM corrections, the exponen- 
tials are expected - from analytic arguments - to be modified by the same 
A— dependent factors as those of ^He. Already for ^^O these factors do not 
deviate much from unity, namely (1 — -j)^^ = y|, (1 — -xzi) = i|- 



4 The Effect of Short-Range Correlations on 



As it has been mentioned above, the evaluation of GMD within the independent- 
particle model is expected to be valid in certain regions of momenta p and Q 
where dynamical correlations do not play a significant role. The next step is 
to consider dynamical, short-range correlations on the GMD. In this section 
we use state-independent central (Jastrow) correlation functions to introduce 



4.1 Low-order approximation for GMD- Application to 



4.1.1 Low-order approximation for GMD 

Our approach is based on the Jastrow formalism and employs the low-order 
approximation (LOA) of Refs. [5H [23, [3S] for the 2DM. Performing the 
spin-isospin summation {v = 4) in Eq.(14) of Ref. [36], we obtain for the 
correlated 2DM p2(n_, "^2; ^i', ^2') the following expression 



GMD 



the SRC. 



4He 



LOA 



ri,r2;ri',r2' 



[1 + g{ri2, rv2')]p2{ri, f2; fr,r2') 
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- jy ^ I / fi'(r34,r34){P2h(^2,r4;f3)p^(fi,f3;fr,f2') 

+ Piin, r3)[pi{r2, ^2')P2h(^3, n; 

- i^-V?(^"^2,rl)p°(f3,r4;ri.,f20]}dV3dV4 , (22) 

where rij = \ fi—fj\, g{r, r') = f{r)f{r') — l and /(r) is the Jastrow correlation 
function, which has to obey the conditions /(O) < 1 and /(r) —>■ 1 for r — »• 
oo. In Eq. fl22|) the uncorrelated IDM p\{f,f''), the 2DM P2(^^5 ''^S ^i'^ "^2') 
and the half-diagonal 2DM p2h(^^5 ''^2; ^i') are calculated in the harmonic- 
oscillator model. The correlated GMD n^'^^{p,Q), is then calculated by 
Fourier transforming according to Eq. (jlj) . The correlated momentum distri- 
bution n^°'^(p) is calculated likewise in LOA by performing the spin-isospin 
summation in expression (13) of Ref. [3B] for the one-body density matrix 
and by Fourier transforming with respect to fl — fi'. Since the LOA preserves 
the normalization of the density matrices, the GMD calculated in the LOA 
obeys the sequential relation ([5]) if also n(p) is calculated within LOA. 

4.1.2 Application to ''He 

The calculation has been carried out for the GMD in the nucleus ^He. The 
momentum distribution n^{p) and the GMD n^{p,Q;b) in the harmonic- 
oscillator model are given by Eqs. fl2Tl) and (!20|) respectively. The HO 
parameter b = 1.382 fm reproduces the experimental value of the charge 
root mean square radius (rms) of ^He, ('"ch.exp)^^^ = 1-67 fm [16| in this 
model if corrections due to the center-of-mass motion and finite nucleon size 
are taken into account. The evaluation of the momentum distribution and 
of GMD in LOA has been carried out using a single- Gaussian correlation 
function, /(r) = 1 — ce~^ . The expression of the GMD so obtained is 

n'^''\p,Q) = n%p,Q-b) + -^[n^'\p,Q-b,c,y) + n^ , (23) 

where y = and 



TV HV Q- b C V) = e 2 l+2y p 4 (l+2j/){l+4j,) 

{P,^,o,c,y) [(i + 2y)(l + %)]3/2^ 



;i + 3y)3/2 



_1 l+2y/ 2^2, 2^2-) _ll+4j/r)2^2 I 2y ffQf,2 
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n^^\p,Q;b,c,y) 



(1 + 3y)3/2 
10c 



e 2 i+3a 



21/ 

e i+3a 



«)-Q62 



+ 



(l + 2?/)3/2 (l + 4?/)3/2 
[(l + 2y)(l + %)]3/2* 



l+2y g (1+2h)(1+4j/) 4 



2c 



+ 



;i + 3y)3/2 

2c 

;i + 32/)3/2 
2c2 

1 + %)3/2 

4c 

;i + 2y)3/2 



_ll+2y. 2^2, 2^2', _1 1 + y r)2j^2 V_g.QI,2 



_1 l+2y(. 2j2, 2^2-) _1 1+;/ n2b2 2/ ^.Qfy2 
g 2 l+3i/'^P ° ° ''g 4 1+3;/^ ° gl+3a"^ 



g-|(p^fe^W62)g-ii±fQ262 



e 2 



i(p262+^2^2) l^Q2fe2 



(24) 



We realize that the n{p^Q) evaluated in LOA fulfils the properties ([5]), ([6]). 
The corresponding expression for the momentum distribution n^'^^{p) of ^He 

[IHlEg] is 



n 



LOA/ 



+ 



+ 



n°{p;b) + 



463 



vr 



3/2 



463 
6c 



c e 



21/ 



2ce" 



D262l±2y 



i+3a 



[(l + 2?/)(l + 4|/)]3/2 

3c2 



3y)3/2 



:i + 2?/)3/2 (1 + 4^)3/2^ 

, _p2{,2l+2y 

4ce ^ i+si' 



[(l + 2y)(l + 4|/)]3/2 (l + 32/)3/2 



(25) 



The term n^'^\ Eqs.f l23p . fl24|) . does not contribute if only the first two terms 
of the LOA approximation are used. This latter case has been considered in 
Ref. HH. 

The inclusion of CM corrections to the correlated GMD is a rather tedious 
task. Only recently, CM and short-range correlations (within the LOA) have 
been considered simultaneously, for one-body quantities, namely the density, 
form factor and momentum distribution of ''He ^45j. In particular, the cor 
responding expression for n{p) if the correlation function /(r) = 1 — ce 
is used for c = 1 is 



n 



LOA+CM 



463 

3/2 



TT 



1 + 



6 



;i + 2y)3/2 (i + 4?/)3/2; 

12 



'_^3/2 



e 3^ 



2^2 



(26) 



24 _i£!ti 48 4(i+2a) 0,2 

_|_ p (3+8:/) e (3+10m)^ " 

[(1 +4y)(3 + 8y)]3/2'' (3 + 10?/)3/2 

For the purposes of the present work we will consider an approximate 
scheme to estimate the CM effects on the correlated GMD of ^He, based on 
the scaling that the CM corrections seem to introduce to the momentum vari- 
ables, see Sec. 3. We scale the variables p and iv by {^ziY^'^ = (4/3)^/^ and 
the variable Q by {^Y^^ = (2/3)^/2^ so that the new GMD obeys the sym- 
metry property of Eq. (jH]). In addition we multiply the resulting expression 

by (i3t) = (4/3)^/^ to ensure the relation between the GMD for Q = 
and the properly normalized momentum distribution Eq. ([5]). The latter is 
obtained from the correlated momentum distribution, as calculated within 
the LOA, by applying the same scaling to the variable p and multiplying 

the resulting expression by [^zi) ■ (We will call the above approximate 
scheme 'LOA+CM'). 

Numerical results for GMD have been obtained with two parametriza- 
tions of f{r). First, the Gaussian correlation function (IG) which has been 
used in evaluating several quantities of ^He [351 EHl EH EH HO, SH [33] has 
been considered. The values of its parameters c and (3 are 0.76 and 0.83 fm 
respectively and the wound parameter k { k = /(I — f{ri2)Yp2{ri2)(Pri2, 
P2{fi2) is the relative pair density distribution function calculated in the HO 
model and normalized to unity) equals 0.018. Second, the Gaussian correla- 
tion function ((72) which has been used for the calculation of GMD of infinite 
nuclear matter with the FHNC/0 approach [2], [18] has been used aiming on 
one hand to compare the LOA results in ^He with the FHNC/0 results m in- 
finite nuclear matter obtained with the same correlation function, and on the 
other to test the sensitivity of the results on the strength of the correlation 
function used. In addition, g2 has been used to estimate the CM corrections 
on n{p), as discussed later. The involved parameters c and (3 equal to 1 and 
0.6765 fm respectively. The wound parameter for g2 equals 0.016. The two 
correlation functions are plotted in Fig. 1. 

Some numerical results for n{p,Q) of ^He are presented in Figs. 2-6. In 
Fig. 2, the GMD is plotted for Q parallel to p, Q = Qpp as a function of 
Qp for p = 0, 1, 2 and 3 fm~^. The correlated GMD calculated with the use 
of LOA and of the correlation functions IG and g2, as given by Eq. (12^ . 
is plotted by continuous and dashed lines respectively. The GMD in the 
harmonic- oscillator model, Eq. fl20l) . is plotted by dotted lines. It seems 
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that SRC within LOA contribute a negative part to n{p^ Q) for positive Qp 
mainly at high values of p. As expected, deviations from the HO picture 
are larger for high values of p and (or) Qp. We should recall that in infinite 
nuclear matter it was found that SRC are mainly significant for \p — Q\ > kp 
and/or p > kp |2]. Results obtained with correlation functions IG and g2 are 
similar, but one should bear in mind that they imply about the same strength 
of SRC (as judged by the corresponding values of the wound parameter). 

Initially the effect of the CM motion in "^He have been considered on the 
generalized momentum distribution per pair for Q = (which equals the 
momentum distribution per particle (Eq. (5)). In Fig. 3 the momentum 
distribution n{p) per particle has been plotted in the single HO model (Eq. 
(21)) and in the HO including CM effects (HO+CM, Eq. (19)) as well as in 
LOA including SRC via correlation function g2 (LOA, Eq. (25)). CM effects 
have been evaluated exactly (LOA+CM, Eq. (26)) and using approximation 
'LOA+CM'. We realize that in ''He in the LOA+CM evaluation of n{p) a 
shrinking is observed relative to its values in LOA which affects the low- 
medium range of p, while the high momentum tail is not affected. A similar 
shrinking occurs in the HO evaluation. In addition one notices that our 
approximation 'LOA+CM' is quite satisfactory for estimating the CM effects 
on n{j)) and can be used also in the case of n{p,Q). In Figs 4 and 5 the 
GMD is plotted for Q parallel p and for Q perpendicular to p for this case, 
'LOA+CM' approximation, as well as for LOA approximation (LOA, Eq. 
(23)) with correlation function g2. For comparison, plots are shown for the 
simple harmonic oscillator (HO, Eq. (20)) and with CM effects (HO+CM, 
Eq.(17)). We realize that CM effects modify the range of values of n{p,Q) 
and produce a similar change in LOA and HO evaluations. 

In Fig. 6, a comparison is made for the GMD per particle for Q par- 
allel to p in the case of ^He (continuous line) and infinite nuclear matter 
(diamond chain) calculated within the LOA and FHNC/0 approximation [2] 
respectively. The Fermi wave number kp of infinite nuclear matter has been 
taken equal to 1.3925 fm~^ (density = 0.182 fm~^) and in both systems 
the correlation function g2 has been used. Results are shown for p equal 
to and 2k^^. We observe qualitative and in some cases even quantita- 
tive agreement in the range where SRC are expected to dominate {p > kp 
and/or \p — Q\ > kp). This implies that the convergence of LOA is satisfac- 
tory. The disappearance of discontinuities in the behavior of GMD of ''He at 
Qp = \p ~ kp\ 01 Qp = \p + kp\ characterizes the transition from infinite to 



14 



finite Fermi systems. 

Regarding the effect of other than central correlations an estimate can 
be drawn for the GMD per pair for Q = (which equals n{p)/A, Eq. (5)). 
In Fig. 3 the results of calculations of n{p)/A with variational Monte Carlo 
method [19] (dashed line) are plotted along with our results. It seems that 
the deviations between our results (in LOA approximation) and those of ref. 
[49] originate mostly from such correlations. 

4.2 Approximation for heavier nuclei using LOA for 
^He. Application to ^^O 

4.2.1 Approximation for heavier nuclei 

In principle, one could evaluate the effect of SRC on the GMD of all £-closed 
shell Z = N nuclei within the LOA approximation for the 2DM (Eq. (1221) ) 
using the results for the one and two-body density matrices in the harmonic 
oscillator model [3] and Fourier transforming according to Eq. (jl]). The 
resulting expressions even in the case of the nucleus ^^O are rather long 
and complicated. We propose an approximation for calculating the GMD 
for not too heavy nuclei which makes use of the GMD of ^He calculated in 
LOA. It has been derived from a corresponding approximation that seems 
to be valid for the momentum distribution n(p). Microscopic calculations of 
n{p) indicate that for large values of p [p > kp) the momentum distribution 
per nucleon is mainly dominated by the SRC and is almost independent 
of y4, whereas for small values of p it is satisfactorily described within the 
independent particle model [HI [50l [51]. Recently, it has also been shown 
experimentally that the momentum distribution at high momenta has the 
same shape for all nuclei differing only by a scale factor [5]. Therefore, 
the following approximation for n{p) seems reasonable for a nucleus of mass 
number A, if we use the momentum distribution per particle of the nucleus 
^He {A = 4) 

n^°"-\p-, A) = n\p; A) + ^An{p; 4) , (27) 

where n^{p;A) and n'^°"~'^{p; A) are the momentum distribution in the in- 
dependent-particle model and including SRC respectively and An{p; 4) = 
n'^°^^{p;4:) — 'nP{p;4)). Using the sequential relation, Eq. ([5]), and the corre- 
lated GMD of the nucleus ^He, we can derive the corresponding approxima- 
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tion for GMD of the nucleus with mass number A 



Q- A) = n\p, Q- A) + ^^^^Anip, Q; 4) , (28) 

where n^{p,Q; A) and n'^°'^'^~'^{p,Q] A) are the GMD in the independent- 
particle model and including SRC respectively and An{p, Q; 4) = n'^°"(p, Q; 4)- 
n°(p, Q; 4). We will use LOA to evaluate n'^^'^ip, Q; 4) (Eq. (28)) and we will 
call the corresponding expression of n'^°'^'^ [p, Q; A) as LOA-4. The approxi- 
mation fl28p is valid for not large values of A, as it does not produce a correct 
asymptotic behaviour for A ^ oo. 



4.2.2 Application to ^^O 

In this paper, we have applied approximation (1251) to calculate the effect of 
SRC on the GMD of the nucleus ^^O. We have considered the special case 
that p and Q are parallel. Results for the GMD of this nucleus within the 
harmonic-oscillator model, n°(p, (5;16), have been presented and discussed 
in Ref. |3]. As for A?7,(p, Q;4) we have used our results for the nucleus 
^He within the LOA for the correlation function IG, presented in section 
4.1.2. The harmonic oscillator parameter b for both nuclei ^He and ^^O has 
been determined in such a way as to reproduce the experimental value of 
the charge r.m.s. radius ((?"ch,exp)^^^ — 1-67 fm and 2.737 fm respectively 
|46j). We found the values 1.2195 fm and 1.7825 fm respectively within the 
harmonic-oscillator model. 

Some of our results are plotted in Figs. 7 and 8. Fig. 7 illustrates the 
variation of n{p, Qp) as a function of Qp for p = 0, 1.1, 2 and 3 fm"^^. Among 
the values of p considered, we have included the value 1.1 fm~^, which cor- 
responds to the Fermi wave number fc^ in ^^O P, The results obtained 
with the approximation (l28l) and within the HO model are displayed. We 
realize, as expected, the effect of SRC mainly for p > kp- Since for the 
description of SRC in the GMD of ^^O use is made of the above evaluation 
of GMD of ^He, the conclusions drawn for their effect on the behavior of 
GMD of ^^O are similar. Fig. 8 provides an estimate of the quality of the 
approximation fl28l) . of the magnitude of the omitted higher-order terms and 
of the contribution of other than central correlations. The GMD per pair for 
Q = calculated within fl28l) . as described above (which equals the momen- 
tum distribution per particle n{p)/A according to Eq. ([5])) is compared with 
the n{p)/A calculated within variational Monte Carlo method [4S] and the 
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harmonic-oscillator model. Judging from the quite good agreement between 
LOA and FHNC results for the momentum distribution found in ref. [12] 
with the use of the same correlation function and the results presented in ref. 
[43j it seems that the major part of the deviation between our results and 
those of ref. [19] stems from other than central correlations. As mentioned 
in Sec. 3 the CM corrections on the GMD of ^^O estimated by using suitable 
scaling for p, w and Q, or a more exact method, are expected to be small. 

5 Summary and Conclusions 

In summary, the study of the generalized momentum distribution n{p, Q) of 
finite, Z = N , ^-closed shell nuclei in their ground state that has been started 
in Ref. [3] within the context of the independent-particle shell model using 
harmonic-oscillator wave functions, was continued by including Jastrow-type 
correlations for investigating the effect of short-range correlations, which is 
expected to be important in certain regions of momenta p and Q. First, the 
low-order approximation of Ref. [SB] has been used and the n(p, Q) of "^He has 
been evaluated using single- Gaussian correlation functions (parametrizations 
IG and g2). Significant deviations from the independent-particle picture 
were found for rather large values of p and (or) Q [p > kp, \p — Q\ > kp). 
The convergence of LOA was explored by comparing with the results of 
the FHNC/0 calculation of n{p,Q) of infinite nuclear matter [2] in which 
the same correlation function g2 was used and it was found satisfactory in 
most cases. In addition, the effects of CM motion has been estimated and 
found not negligible and the role of correlations other than central has been 
brought up. Next, an approximation scheme for the evaluation of n{p, Q) 
of heavier nuclei was proposed (Eq. fl28l) ) that includes the effect of short- 
range correlations by means of the above evaluated GMD of ^He. Numerical 
results have been derived for ^^O and the quality of the approximation was 
discussed. 

Further investigation of n{p, Q) of finite nuclei should consider the exact 
evaluation of the corrections due to the Center-of-Mass motion which are 
quite significant in the case of light nuclei as we have realized in the ap- 
proximate treatment of ''He. One should start along the lines of references 
[HmS] for the case of N = Z, 1-closed nuclei. Starting from these nuclei, one 

— * 

must also consider to include in the calculations of n{p, Q) state-dependent 
correlations. Another interesting direction for future work is the determina- 
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tion of other Fourier transforms of the two-body density matrix, for example 
n{p,k,Q) [521 The above evaluation of n{'p,Q) is a first step towards 
more realistic calculations which include also other than Jastrow correlations. 
The quantity n(p, Q) is mainly useful for the study of final-state interactions 
of struck nucleons as they propagate inside the nuclear medium in various 
scattering processes and for the understanding of elementary excitations of 
nuclei. 
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Figure 1: The single-Gaussian correlation functions with c = 0.76, (3 = 0.83 
fm (IG) (continuous line) and with c = 1, (3 = 0.6765 fm (5(2) (dashed line) 
(see 4.1.2). 
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Figure 2: Generalized momentum distribution of "'He for Q parallel to p, 
n{p, Qp) as a function of Qp for p = 0, 1, 2, 3 fm^^ including SRC with LOA, 
Eq.(23) using correlation functions IG (continuous line) and g2 (dashed line) 
and in the harmonic-oscillator model, Eq. (20) (dotted line), (see 4.1.2) 
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Figure 3: The momentum distribution per particle (or the generahzed mo- 
mentum distribution per pair for Q = (Eq.(5)) as a function of p in the 
case of ^He in the simple harmonic oscillator model (HO, Eq.(21)) and with 
CM effects (HO+CM, Eq. (19)) as well including SRC using LOA and cor- 
relation function g2 (LOA, Eq.(25)) and considering CM effects (LOA+CM, 
Eq.(26)) and approximation 'LOA+CM' (see sec. 4.1.2). In addition, the 
results of the variational Monte Carlo calculation VMC [49] are plotted. A 
logarithmic scale is used. 
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Figure 4: Generalized momentum distribution of He for Q parallel to p, 
n{p, Qp) as a function of Qp for p = 0, 1, 2, 3 fm~^ in the simple harmonic 
oscillator (HO, eq.(20)) and with CM effects (HO+CM, eq. (17)) as well 
including SRC using LOA and correlation function g2 (LOA, eq.(23)) and 
considering CM effects using approximation 'LOA+CM' (see 4.1.2). 
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Figure 6: Generalized momentum distribution per particle for Q parallel to 
p, n{p, Qp)/A as a function of Qp for ^He calculated by including SRC using 
LOA (continuous line, Eq.(23)) and for infinite nuclear matter with Fermi 
wave number kp^^ = 1.3915 fm~^ calculated in Fermi-liypernetted chain 
approximation [2] (diamond chain). In both systems the correlation function 
g2 is used and results are presented for p = and 2 kp'^ . 
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Figure 7: Generalized momentum distribution of ^^O for Q parallel to p, 
n(p, Qp) as a function of Qp for p = 0, 1.1, 2 and 3 fm~^ including SRC using 
the approximation of Eg. (1281) and LOA (LOA-4) and correlation function IG 
(continuous line) and in the harmonic-oscillator model (dotted line). 
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Figure 8: Generalized momentum distribution per pair for Q = (or mo- 
mentum distribution per particle (Eq. (5)) as a function of p, in the case of 
^^O, including SRC using the approximation of Eq. fl25]) and LOA (LOA-4) 
and correlation function IG (continuous line) along with the results for the 
momentum distribution per particle of a variational Monte Carlo calculation 
[39] (dashed line, VMC) and of the harmonic-oscillator model (dotted line). 
A logarithmic scale is used. 
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